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The process pp — + ttbb + X represents a very important background reaction to searches at the LHC, in 
particular to ttH production where the Higgs boson decays into a bb pair. A successful analysis of ttH at the 
LHC requires the knowledge of direct ttbb production at NLO in QCD. We take the first step in this direction 
upon calculating the NLO QCD corrections to the subprocess initiated by qq annihilation. 



1. Introduction 

The search for new particles will be the primary 
task of the LHC. The discovery of new particles 
first of all requires to establish excess of events 
over background. The situation at the LHC 
is particularly complicated by the fact that for 
many expected signals the corresponding back- 
ground cannot entirely be determined from data, 
but has to be assessed upon combining measure- 
ments in signal-free regions with theory-driven 
extrapolations. To this end, a precise prediction 
for the background is necessary, in particular in- 
cluding NLO corrections in QCD. Since many of 
these background processes involve three, four, or 
even more particles in the final state, this kind of 
background control requires NLO calculations at 
the technical frontier. This problem lead to the 
creation of an "experimenters' wishlist for NLO 
calculations" at the Les Houches workshop 2005 
[1], updated in 2007 [2], which triggered great 
theoretical progress in recent years (see ReL [2] 
and references therein). Meanwhile the listed pro- 
cesses involving at most five particles in loops 
have been completed in NLO QCD, including the 
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production of WW-|-jet [3], weak-boson pairs plus 
two jets via vector-boson fusion [4], and triple 
weak-boson production [5,6]. However, none of 
the true 2^4 processes has yet been addressed 
at NLOo Among those processes, pp -^ ttbb -I- X 
has top priority. 

This process represents a very important back- 
ground to ttH production where the Higgs boson 
decays into a bb pair. While early studies of ttH 
production at ATLAS [8] and CMS [9] suggested 
even discovery potential of this process for a light 
Higgs boson, more recent analyses [10, 11] with 
more realistic background assessments show that 
the signal significance is jeopardized if the back- 
ground from ttbb and tt -|- jets final states is not 
controlled very well. This is a clear call for im- 
proved signal and background studies based on 
NLO predictions to these complicated processes. 
For the ttH signal [12] and the tt -I- Ijct back- 
ground [13] at the LHC such predictions have 
been accomplished in recent years. 

The dominant mechanism to produce ttbb fi- 
nal states in hadronic collisions is pure QCD. 
In leading order (LO) quark-antiquark (qq) and 
gluon-gluon (gg) initial states contribute, where 
the latter strongly dominate at the LHC because 
of the high gluon flux. Being of order a^ the cor- 
responding cross sections arc affected by a very 
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the following way, 



Figure 1. Sample diagrams contributing to qq 
ttbb in LO QCD. 





Figure 2. Pentagon and hexagon sample graphs 
contributing to qq — > ttbb at one loop in QCD. 



large scale uncertainty, which amounts to a fac- 
tor two or more. Technically the qq channel is 
simpler to deal with — though still demanding — 
and thus represents a natural first step towards 
a full treatment of pp — > ttbb + X at NLO. In 
these proceedings we briefly report on this first 
step, which was accomplished in Ref. [14]. 

In LO QCD seven different Feynman dia- 
grams (Fig. [T]) contribute to the partonic process 
qq — > ttbb. The virtual QCD corrections com- 
prise about 200 one- loop diagrams, which include 
8 hexagons and 24 pentagons (Fig. [5]). The real 
corrections, qq -^ ttbbg, are described by 64 di- 
agrams, which result from the LO graphs upon 
adding an external gluon in all possible ways. In 
the following we briefly describe the calculation 
of the virtual and real corrections and present 
numerical results for the LHC. The calculation 
has been worked out twice and independently, re- 
sulting in two completely independent computer 
codes, which we denote as version 1 and 2. 

2. Virtual corrections 

The general strategy for the evaluation of the 
one-loop corrections is based on the reduction of 
the amplitude Al^'"-' of each (sub)diagram F in 



M(r)^C(i^) 






(1) 



where the colour structure C^^' present in the 
(sub)diagram is factorized from the remaining 
colour- independent part. While loop diagrams 
that involve a quartic gluon vertex give rise to 
3 colour- factorized amplitudes A4^^\ in all other 
diagrams — which constitute the large majority — 
the colour structure factorizes completely. This 
implies that the computation time for individual 
loop diagrams does not scale with the total num- 
ber of independent colour structures, which is 6 
for qq -^ ttbb. 

The colour- free part of A^ f'"' is written as a lin- 
ear combination of so-called standard matrix ele- 
ments (SMEs) Mm, which consist of Dirac chains 
containing the polarization information. Since 
the computing time scales with the number of 
SMEs, it is important to reduce the set of SMEs 
as much as possible. To this end, we employ an 
algebraic procedure based on four-dimensional re- 
lations that are derived from Chisholm's identity 
and are applied after separation of UV and IR 
singularities. For massless external fermions this 
four-dimensional reduction has been described in 
detail in Sections 3.1 and 3.3 of Ref. [15]; here 
we had to generalize this approach to one mas- 
sive and two massless spinor chains (b quarks are 
treated as massless particles). In this way some 
thousand different spinor chains are reduced to 
about 150 SMEs Mm- 

The most time-consuming components of the 
numerical calculation are the scalar form fac- 
tors J-"m , which are linear combinations of the 
Lorentz-invariant coefficients of A^-point tensor 
loop integrals with rank R < 3 and degree N < 6, 

:fP = V V ic^Il ..T^ ... (2) 
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The evaluation of the tensor coefficients T^ ■ 
follows the strategy of Refs. [16, 17] that was al- 
ready successfully used to compute the NLO elec- 
troweak corrections to e+e" -^ 4 fermions [15,18]. 
In this approach the analytic expressions are not 
reduced to master integrals. In contrast, the ten- 
sor integrals are evaluated by means of algorithms 



that perform a recursive reduction to master in- 
tegrals in numerical form. This avoids huge an- 
alytic expressions and permits to adapt the re- 
duction strategy to the specific numerical prob- 
lems that appear in different phase-space regions. 
The scalar master integrals are evaluated using 
the methods and results of Refs. [19,20]. UV di- 
vergences are regularized dimensionally in both 
evaluations, but IR divergences are treated in 
different ways as described below. Tensor and 
scalar 6-/5-point functions are directly expressed 
in terms of 5-/4-point integrals [16, 17]. Ten- 
sor 4-point and 3-point integrals are reduced to 
scalar integrals with the Passarino-Veltman algo- 
rithm [21] as long as no small Gram determinant 
appears in the reduction. If small Gram deter- 
minants occur, two alternative schemes are ap- 
plied [17]. One method makes use of expansions 
of the tensor coefficients about the limit of vanish- 
ing Gram determinants and possibly other kine- 
matical determinants. In the second (alternative) 
method we evaluate a specific tensor coefficient, 
the integrand of which is logarithmic in Feynman 
parametrization, by numerical integration. Then 
the remaining coefficients as well as the standard 
scalar integral are algebraically derived from this 
coefficient. 

Version 1 of the virtual corrections starts with 
the generation of Feynman diagrams using Fey- 
nArts 1.0 [22]. Their algebraic reduction is 
completely performed with in- house Mathemat- 
ICA routines. In detail, D-dimensional identities 
(Dirac algebra, Dirac equation) are used until UV 
divergences cancel against counterternis. The IR 
(soft and coUinear) divergences are regularized 
dimensionally and expressed in terms of scalar 
2- and 3-point integrals as described in Ref. [23], 
keeping the full dependence on _D. It turns out 
that no D-dependent prefactors occur in front of 
the IR-singular i?o and Co integrals. In Ref. [14] 
we show that this result of our specific calculation 
is not accidental, but generalizes to arbitrary pro- 
cesses at NLO in the Feynman gauge. Having 
cancelled UV divergences against counterterms 
and controlled the I?-dimensional issues concern- 
ing IR singularities, the amplitude is further sim- 
plified by reducing the SMEs in four space-time 
dimensions. 



Version 2 of the virtual corrections employs 
FeynArts 3.2 [24] for generating and FORM- 
Calc 5.2 [25] for preprocessing the amplitudes. 
The first part of the calculation is performed in D 
dimensions. In particular, the so-called rational 
terms resulting from the UV divergences of ten- 
sor loop coefficients are automatically extracted 
by FormCalc. Since the IR divergences that 
appear in the qq channel are of abelian nature, 
we exploit the fact that they can be regularized 
as in QED by means of fermion and gauge-boson 
(gluon) masses, mq and Trig. These masses are 
treated as infinitesimal quantities (with m^ <C 
mq) both in the algebraic expressions and in the 
numerical routines that evaluate the tensor in- 
tegrals, i.e. only the logarithmic dependence on 
these mass parameters is retained. Correspond- 
ing IR singularities associated with real emission 
have been obtained from Ref. [26] by means of an 
appropriate change of regularization scheme. 

3. Real corrections 

In both evaluations the singularities from soft 
or coUinear gluon emission are isolated via dipole 
subtraction [26 28] for NLO QCD calculations 
using the formulation [26] for massive quarks. Af- 
ter combining virtual and real corrections, sin- 
gularities connected to coUinear configurations in 
the final state cancel for coUinear-safe observables 
automatically after applying a jet algorithm. Sin- 
gularities connected to coUinear initial-state split- 
tings are removed via factorization by PDF redefi- 
nitions. The phase-space integration is performed 
with multichannel Monte Carlo generators [29] 
and adaptive weight optimization similar to the 
one implemented in RacoonWW [30]. 

In version 1 of the real corrections the ma- 
trix elements have been calculated using the 
Weyl-van-der-Waerden spinor technique in the 
formulation of Ref. [31]. Soft and coUinear sin- 
gularities are regularized using dimensional reg- 
ularization. The phase-space integration, im- 
plemented in C-|--f , is based on RacoonWW, 
but the phase-space mappings are build up in a 
more generic way very similar to the approach of 

LUSIFER [32]. 

In version 2 of the real corrections the ma- 



(7L0 [fb] PP ^ "^*^ + ^ 











hh "T yy 






000 




^s^^ 


100 




^^^*^ 


10 


■ nit = 


172.6 GcV ^V^^^H 




«)/2 


< /t < 2^0 



20 40 60 80 100 120140 160 180 200 

™bb,cut[GeV] 



™^ - 1 [%] pp(gg) ^ ttbb + X 



10 



1 1 1 1 r— 

slicing t— X— -! 
subtraction 

mt = 172.6 GeV 

<5c = 10-« 



_l I I I L_ 



-5 -4 -3 -2 -1 
logio k 



Figure 3. Complete LO cross section (gg + qq) 
and qq contribution as a function of mj^g ^^j. 



Figure 4. NLO corrections to Y'Vill) ~^ ttbb + X: 
dipole subtraction versus phase-space slicing. 



trix elements have been generated with Mad- 
GRAPH [33]. As in the virtual corrections, soft 
singularities are regularized by an infinitesimal 
gluon mass and coUinear singularities by small 
quark masses, which appear only in logarithms 
in the endpoint part of the subtraction function. 
The Monte Carlo generator is a further develop- 
ment of the one used in COFFER77 [34] and 
for the calculation of the NLO corrections to 
PP ^^ JJH + X [35]. For the purpose of check- 
ing, we have also implemented two-cut-off slic- 
ing. In this approach the singular parts of the 
phase-space integration, which originate from the 
soft region Eg < Ss\/S/2 and the coUinear regions 
1 — cos(6'gq) < (5c, are treated by means of the 
well-known soft/coUinear approximations. When 
adding all contributions, the 0{ds,c) dependence 
on the technical cuts cancels if the cut-off param- 
eters are chosen small enough. 

4. Numerical results 

We present predictions for pp — > ttbb -I- X 
at y/s = 14TeV. For the top-quark mass we 
take mt = 172.6 GeV, while all other QCD par- 
tons (including b quarks) are treated as mass- 
less. Final states involving coUinear gluons and 
b-quarks are recombined into coUinear-safe jets 



by means of the fcx-algorithm of Rcf. [36] . Specif- 
ically we require t wo b-quark jets with separa- 
tion -^/Af^^ -(- Ay2 > D = 0.8 in the rapidity- 
azimuthal-angle plane. This sets an effective 



lower limit on the invariant mass 



'bb 



of the 



b-quark pair, thereby avoiding coUinear singu- 
larities from g ^ bb splittings. Motivated by 
the search for a ttH{H -^ bb) signal at the 
LHC [10, 11], we impose the following additional 
cuts on the transverse momenta, rapidity and in- 
variant mass of the b-quark jets: px.b > 20 GeV, 
[j/b] < 2.5 and m^g > mbb.cut- We use CTEQ6 
PDFs [37] with A^° = 165 MeV at LO and 
A|*^ = 226 MeV at NLO. For the renormalization 
and factorization scales we use the central value 
/io = m,t -I- m^^b cut/2 and, to estimate the scale 
uncertainty, we vary these scales in the range 
/io/2 < /i < 2/io in a uniform and antipodal way, 
i.e. with /ip{ = ^p = ^ and /xr = /ig/MF = M- 

Figure [3] shows the total LO cross section and 
the contribution induced by qq annihilation as 
a function of rribS cut • The gg channel domi- 
nates over qq annihilation by roughly a factor 
of 17. The scale uncertainty, displayed by the 
bands, amounts to a factor 1.6 and is domi- 
nated by the /iR-dependence owing to the large 
power of as(^R)^ in the LO cross section. Fig- 
ure [4] illustrates the mutual agreement between 
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Figure 5. LO and NLO scale dependence of 
PP('?9) ~^ ttbb + X (antipodal scale variation). 



Figure 6. LO and NLO pp{qq) - 
section as a function of ttij^j, ^.^j. 



ttbb + X cross 



NLO results for the qq channel obtained with 
dipole subtraction and phase-space slicing. Here 
we set rrii^g ^jjj = 0, i.e. /^o = rrit- We ob- 
serve that, within integration errors, the slicing 
results become independent of the soft cut-off Ss 
for 6s ^ 10^"^ and agree with the result in the sub- 
traction method. For the (5c-dependence we find a 
similar behaviour [14]. Figure [5] shows the reduc- 
tion of the scale dependence of the qq contribution 
upon going from LO to NLO. Setting m^^^ ^^^ — 0, 
at /i = /io = TTit we find ctlo = 85.522(26) fb 
and tJNLO = 87.698(56) fb, and we observe that 
(uniform and antipodal) scale variations in the 
range /xo/2 < jjl < 2/io shift the cross section by 
55% in LO and by 17% in NLO. Figure E] dis- 
plays the dependence of the LO and NLO cross 
sections on the invariant-mass cut m-bScut- The 
central scale and the uncertainty bands are as 
in Figure [H The reduction of the scale uncer- 
tainty from about ±50% to ±17% holds true 
for the considered range in Ti^b.cut- While the 
NLO prediction is consistent with the LO uncer- 
tainty band, the shape of the distribution is dis- 
torted by the corrections. For the central scale, 
/iQ — mt + rjibb cut/2i we find an NLO correction 
of +2.5% for small TTibb cut but a correction of 



-10.6% for 



'bb^cut 



= 200 GeV. 



5. Numercial stability and CPU efficiency 

In order to validate the numerical stability of 
the tensor reduction we compared the results 
of the two independent codes finding very good 
agreement both for single phase-space points and 
integrated quantities. The CPU time needed to 
evaluate the virtual corrections for a phase-space 
point (including sums over colours and polariza- 
tions) amounts to about 10 ms on a 3 GHz Intel 
Xeon processor. This demonstrate that the em- 
ployed diagrammatic techniques permit a fast and 
numerically stable evaluation of six-particle pro- 
cesses at the LHC. Based on these encouraging 
results, we expect to be able to extend this calcu- 
lation to the technically more challenging gluon- 
fusion channel. 
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